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Several quantum critical compounds have been argued to have multiple instabilities towards or¬ 
ders with distinct dynamical exponents. We present an analysis of a quantum multicritical point in 
an itinerant magnet with competition between ferro- and antiferromagnetic order, modelled using 
Hertz-Millis theory. We perform a one-loop renormalization group treatment of this action in the 
presence of two dynamical exponents. In two and in three dimensions, when both incipient orders 
are quantum critical, we find that the specific heat, thermal expansion and Griineisen parameter 
obey the same power laws as those expected for a single ferromagnetic quantum critical point. The 
antiferromagnetic correlation length and boundary of the antiferromagnetic ordered phase are sup¬ 
pressed by the dangerously irrelevant interactions with quantum critical ferromagnetic fluctuations. 
We find no difference between a quantum bicritical point and a quantum tetracritical point. Our 
results are compared with experiments on NbFe 2 . 


PACS numbers: 73.43.Nq, 71.10.Hf, 64.40.Kw 

Quantum criticality is characterised by universal diver¬ 
gences of thermodynamic quantities at a continuous zero- 
temperature phase-transition as some non-thermal con¬ 
trol parameter (e.g. pressure, doping or magnetic Held) is 
changed. The much-studied power laws associated with 
the quantum phase transition depend on the universal¬ 
ity class, which in contrast to the classical case depends 
on the dynamics of the order parameter fluctuations in 
imaginary time. The dynamic effects are characterized by 
a dynamical exponent, 2, which, in addition to control¬ 
ling the power-law divergences, defines the boundaries of 
distinct regions in the phase diagram W- 

In recent years quantum critical behaviour has been 
observed in systems that do not seem to conform to the 
standard theoretical picture of fluctuations of an order 
parameter with a unique dynamical exponent. For exam¬ 
ple the specific heat in NbFe2 displays aC~ —TlnT re¬ 
lation as would be expected for three dimensional itiner¬ 
ant ferromagnetic quantum criticality which is described 
with z = 3 . In contrast the resistivity displays p ~ T 3 / 2 
as expected for a dirty three dimensional antiferromag- 
net, usually described with 2 = 2 0E]. 

In this Letter we report results of a study into quan¬ 
tum criticality where an itinerant material is unstable to 
both ferro- and antiferromagnetic order. Materials with 
a multicritical point in the phase diagram [as shown in 
Figs. 1 (a) and 1 (b)] show a competition between two 
distinct types of order. Here we assume that if we had 
control over another tuning parameter we could suppress 
this multicritical point to zero temperature, to form a 
quantum multicritical point [as shown in Fig. 1 (c)]. 
This type of phase diagram has been considered before 
[ 7 ], but in the presence of a symmetry-breaking field, 
which we do not treat here. To model quantum mul¬ 
ticriticality we construct an effective action in terms of 
spin-fluctuations by adapting Hertz-Millis theory mmm- 
This neglects important non-analytic terms which some¬ 
times arise from integrating out the fermionic degrees 


of freedom Jj. Ignoring such terms seems valid slightly 
away from the critical point [ 9 ], Nevertheless their inclu¬ 
sion could stabilize a finite-momentum spin-density wave 
near the ferromagnetic quantum critical point m to gen¬ 
erate the scenario considered here. We extend Millis’ 
calculation [l] to treat the quantum multicritical point 
by following the renormalization group (RG) procedure. 
This enables us to map out the rich phase diagram and 
calculate the leading order critical parts of the specific 
heat, thermal expansion and Griineisen parameter. 

Our main result is that, since we are at or above the 
upper critical dimension of our model (d+ = 2), we can 
essentially treat the fluctuations of each type of order in¬ 
dependently, and the critical part of the free energy and 
therefore its derivatives simply become the sum of the 
contributions associated with each individual order pa¬ 
rameter. The caveat is that whenever the dangerously 
irrelevant interactions affect the RG flow, due to the 
multiple dynamics they can produce novel temperature 
dependences. This allows the dangerously irrelevant in¬ 
teractions to shape the phase diagram. We also discuss 
the resistivity in relation to experiments, though do not 
offer an explicit calculation here. 

Our analysis is formulated in terms of an order param¬ 
eter field 4 >(q,uj n ) which describes the magnetization of 
the system. In a system unstable towards both ferro- 
and antiferromagnetic order, the susceptibility will be 
large near zero momentum and near the antiferromag¬ 
netic wavevector Q. We split the magnetization field 
into two parts to represent these different regimes. We 
denote the small momentum part of the field as 03 (q,u> n ) 
for \q\ < A3 where the subscript 3 refers to the fact that 
this field is usually described by a dynamical exponent 
2 = 3 . To denote the field near momentum Q we use the 
notation 02(<7,w n ) where here q (which we restrict such 
that | q | < A2) measures the deviation from Q. Here the 
subscript 2 refers to the dynamical exponent 2 = 2 usu¬ 
ally associated with antiferromagnetic order. We allow 


2 


T T 




T T 




FIG. 1: Suppressing a multicritical point to zero temperature, 
(a) and (b) show bicritical and tetracritic.al points respec¬ 
tively. By invoking another non-thermal control parameter g , 
these multicritical points can be tuned to zero temperature, 
as shown in (c). (d) shows the expected phase diagram when 
two quantum critical points arise for increasing g. 

for <^3 and </>2 to have 713 and 712 components respectively. 

The action we use to describe a quantum multicritical 
point is the sum of the actions of a ferromagnetic and an 
antiferromagnetic quantum critical point (QCP), 

S [03, fa] = Xi 1 w„) 

i— 3,2 UI n q 

+ j dxdr \u 3 (j)\{x 1 t) + ui<j>\{x, r) 

+ « 32 <Pl( x,t)4%(x,t)\ , ( 1 ) 

where the two QCPs are coupled together by a mode¬ 
mode coupling term 2232 . The bare inverse spin suscepti¬ 
bilities are given by 

= 6i + q 2 + ( 2 ) 

g~* - 

where 23 is the dynamical exponent associated with ferro¬ 
magnetic order, 2 = 3 , and 22 is the dynamical exponent 
associated with antiferromagnetic order, 2 = 2 . We have 
added ‘kinetic coefficients’ 773 and 772 which allow us to 
renormalize in the imaginary time direction [mug. The 
classical analogue of the action describes a multicritical 
point in the <53-52 plane. The model shows bicriticality 
if rt | 2 > 4723722 and tetracriticality if this is inequality is 
reversed mm- 

In order to map out the phase diagram and predict 
thermodynamic quantities, we perform an RG analysis by 
simultaneously integrating out the fa modes in a small 
shell with momenta between A3 jb and A3, and the fa 


modes with momenta between A2/6 and A2. We then 
rescale such that the original cut-offs are restored and 
calculate how the other parameters in the model must 
rescale. The presence of multiple dynamical exponents 
means there is no unique way to rescale frequency. We 
choose to rescale it as b z where 2 is a fictitious dynamical 
exponent which we leave unspecific (see Ref. . TT and 
m)- This enables renormalization but will drop out of 
our calculations so that no physical properties depend on 
it. The RG equations can either be derived directly from 
the action or by calculating a physical property (such as 
the free energy) and ensuring it does not change under 
renormalization. We have done both, and find that the 
one-loop RG equations for the tuning parameters and 
interactions are 

^ = 25,(6) + 22,(6)4(n, + 2)/f> (5,(6),71(6)) 

+ w,( 6 ) 472 ,/, (2) ( 5 ,( 6 ), 7 i( 6 )), ( 3 a) 

= [ 4 - (d + z i )]u i (b)-A{n i +8)f^\d i (b),Ti(b))u^ 
- 4 ti,/,^( 5 ,( 6 ), 7 i( 6 ))w,( 6 )w,( 6 ), ( 3 b) 

|^ = [4 — (<(+*)] #,(<>) 

- 4 K'+ 2 )/j 4) ( <5 i(6)> T j {. b ))u j ( b )w i (b), ( 3 c) 

i= 3,2 

where i is either 3 or 2 and 1 is correspondingly either 2 
or 3 . We have defined 22, = 72 ,/t?, and 76, = 7232/77,, as we 
find that the only way the interactions enter the RG equa¬ 
tions is in these combinations. For the rest of this Let¬ 
ter, when a parameter is written without explicit scale- 
dependence we are referring to the bare, unrenormalized 
value. We find that the temperature scales as T( 6 ) = b z T 
and the kinetic coefficients are 77,(6) = b Zi ~ z r]i. The RG 
equations can be written in terms of two temperature 
fields 7 ,( 6 ) = ?7,(6)T(6) = 7 ) 6 2i which represent the ef- 
fective temperature felt by the fa modes. The f- and 
,/ 44 ' ) functions are the one-loop integrals shown in Fig. 
2 , which arise in the RG procedure due to interactions 
with modes above the cut-off. They are exactly the same 
functions which appear in Hertz-Millis theory, defined in 
Ref. [ 1 ] . The scaling of the interaction terms can be cal¬ 
culated from Eqs. ( 3 b) and ( 3 c), which can be analysed 
under the usual approximation that the functions 
are constant [ 4 ]. In d = 3 we find all interaction terms 
decay as some power of 6 , whereas in d = 2 , 722 and 
762 decay logarithmically. We find that at large values 
of 6, 222(6) ~ (ln6) _1 just as in the 7 = 2 = 2 Hertz- 
Millis case, and 762(6) ~ ^ n 7,^ _ ( n 2+ 2 )/( n 2+8)^ w j 1 j c j 1 a 
very slow decay unique to the multicritical case. Because 
the free energy can be written as a power series in these 
interaction terms, we conclude that the upper critical di¬ 
mension is 7 / = 2 , just as for an antiferromagnetic QCP. 
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FIG. 2 : One-loop diagrams that contribute to the RG flow 
and Eqs. ( 3 a), ( 3 b) and ( 3 b). Lines with free (connected) 
ends represent inodes below (above) the cut-off. The symbols 
i,j, k = 3 , 2 label which modes the lines represent, (a) shows a 
contribution to the renormalization of Si from interaction with 
4>j modes above the cut-off, and is proportional to /j 2 ' ( 8j,Tj )■ 
(b) shows a contribution to the renormalization of Ui if j = i 
and U32 if j A i, due to interactions with k modes above the 
cut-off, and is proportional to fj, 4> ( 8 k, 71 -). 


In the cases of two and three dimensions considered here, 
we are at or above the upper critical dimension and so 
are controlled by the Gaussian fixed point where all in¬ 
teractions flow to zero. 

The distinct regions of the phase diagrams and the 
correlation lengths in each regime can be calculated from 
Eq. 3 (a). In our calculation we use Millis’ approxima¬ 
tion for the integral of the f^(S,T) function, which is 
different in the quantum critical and quantum disordered 
regimes g], defined by R <C T 2 ^ z and R T 2 / 2 for a 
single QCP. Here R is the renormalized tuning parameter 
or quasiparticle mass, which may acquire some temper¬ 
ature dependence. In the multicritical case, we conclude 
that both the ferro- and antiferromagnetic modes can in¬ 
dependently be quantum critical or quantum disordered, 
splitting the phase diagram into four regions separated 
by the two lines Ri ~ T ^ Zi . The solution of Eq. 3 (a) 
at large values of b yields a tuning parameter becomes 
Si(b) = b 2 Ri where Ri is renormalized by interactions 
with both modes independently. This can be related to 
the correlation length of the corresponding order param¬ 
eter by Ri = £~ 2 . We denote the zero temperature part 
of Ri by ri, and it is this which tunes to the QCP at 

n = 0 . 

It is the dangerously irrelevant interaction terms which 
give the renormalized tuning parameters their tempera¬ 
ture dependence, which in turn control the correlation 
lengths and the boundaries of the ordered phases. In 
three dimensions the boundaries of the ordered phases 
can be calculated from the lines T(?y) where the corre¬ 
sponding correlation length diverges. While no true or¬ 
der can exist in 2 D, we adopt the usual convention and 
use the point that the Ginzburg criterion of the classical 
theory breaks down to identify the ‘phase boundary’. 

The generic phase diagram for a quantum multicriti¬ 


cal point is shown in Fig. 3 , which is qualitatively the 
same in both two and three dimensions. We find it most 
revealing to interpret the results as the sum of two quan¬ 
tum critical points. In both two and three dimensions, 
we find that to leading order the ferromagnetic z = 3 
QCP is qualitatively unaffected by the antiferromag¬ 
netic z = 2 QCP. However the antiferromagnetic QCP 
is strongly affected by the proximity to a ferromagnetic 
QCP. When the ferromagnetic modes are quantum criti¬ 
cal, the antiferromagnetic correlation length acquires the 
temperature dependence of the ferromagnetic correlation 
£2 2 ~ r2+AT 4 / 3 instead of the usual ~ V2+BT 3 / 2 in 
three dimensions, and ~ r2 + CTln(l/T) instead of 
the usual = r2 + D In (In ( 1 /T)) / In ( 1 /T) in two di¬ 
mensions. This temperature dependence dominates the 
antiferromagnetic correlation length in region I of the 
phase diagram in Fig. 3 . Interactions with quantum 
critical ferromagnetic fluctuations therefore reduce the 
antiferromagnetic correlation length, which in turn sup¬ 
presses the boundary of the antiferromagnetic phase. If 
the ferromagnetic fluctuations are Fermi liquid-like then 
the antiferromagnetic QCP is qualitatively unaffected. 

Since the interactions are irrelevant, thermodynamic 
properties in the disordered region of the phase dia¬ 
gram can be obtained from the Gaussian part of the 
free energy, which can be calculated directly from the 
action. This is just the sum of the contributions from 
both individual QCPs weighted by 77F , Fg(T,5 3 ,5 2 ) = 

^4 3) (S3,r 3 ) + ^ 4 2) (^,r 2 ), where F^faTl) is the 
Hertz-Millis free energy for order with dynamical expo¬ 
nent Zi defined explicitly in Ref [ 3 ]. While the free en¬ 
ergy is described by the Gaussian part, the effect of the 
interactions is seen in the rescaling of the Gaussian pa¬ 
rameters. 

For a single QCP, the specific heat and thermal ex¬ 
pansion behave differently in the quantum critical and 
Fermi liquid regimes, as tabulated in Ref. m ■ In the 
multicritical case, in each of the four distinct regions of 
the paramagnetic phase in Fig. 3 the specific heat and 
thermal expansion are just the sum of the contributions 
from each individual QCP, which we find to leading or¬ 
der are unchanged by the interactions. In both two and 
three dimensions, in region I of the phase diagram in Fig. 
3 , where both ferro- and antiferromagnetic modes would 
expected be to quantum critical, then the strong tem¬ 
perature dependence of the ferromagnetic contribution 
dominates, and the presence of antiferromagnetic quan¬ 
tum criticality is subleading. In the other regions of the 
disordered phase, the observed quantities are the sum 
of the two terms. When the QCPs are separated, suffi¬ 
ciently close to each QCP the effects of the other QCP 
are not measurable and the thermodynamics is as would 
be expected from a single QCP. 

The exception to this is in two dimensions in region I 
of the phase diagram in Fig. 3 . In this case the tempera¬ 
ture dependent renormalization of the antiferromagnetic 
tuning parameter is strong enough to push the antiferro¬ 
magnetic mode out of the quantum critical regime, as the 
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FIG. 3: Generic phase diagram of a quantum multicritical 
point, derived by setting = g + r and V 2 = g — r. The re¬ 
gions in the paramagnetic regime have been identified. I and 
II are quantum critical regions where the ferromagnetic con¬ 
tributions dominate specific heat and thermal expansion. In 
I the ferromagnetic fluctuations control the antiferromagnetic 
correlation length and the boundary of the antiferromagnetic 
ordered phase. In II the antiferromagnetic correlation length 
is dominated by the tuning parameter. Ill is an antiferro¬ 
magnetic quantum critical region, where the antiferromag¬ 
netic correlation length and boundary of the ordered phase 
are dominated by the antiferromagnetic fluctuations, and the 
ferromagnetic correlation length is dominated by the tuning 
parameter. IV is a Fermi liquid. The black dot indicates 
approximately where the antiferromagnetic phase boundary 
undergoes a crossover from a power law associated with an¬ 
tiferromagnetism to one normally associated with ferromag¬ 
netism. The subleading contributions depend on the dimen¬ 
sionality of space, as explained in the main text. 


condition R 2 < T will never be satisfied. This means we 
must use the Fermi liquid approximation (I ?2 > T) in cal¬ 
culations of the antiferromagnetic contribution to phys¬ 
ical quantities, but the correlation length is still domi¬ 
nated by temperature. However in this regime the ferro¬ 
magnetic contributions dominate specific heat and ther¬ 
mal expansion and this effect is subleading. 

We now compare our theory with the existing experi¬ 
mental results in Nbi_ y Fe 2 + y . Near y = 0 this material 
shows both ferro- and antiferromagnetic quantum critical 
points mm There the measured specific heat shows 
a — TlnT dependence consistent with the dominance of 
ferromagnetic fluctuations as we have shown above. The 
thermal exapansion has not been measured but we pre¬ 
dict it to show a T 1 2 / 3 dependence at low temperatures. 
In this work we have not calculated the resistivity be¬ 
cause of the complex interplay we anticipate between hot- 
spot/line scattering of the antiferomagnet [T8] and the 


small angle scattering for the ferromagnetic fluctuations. 
The measured resistivity is Ap ~ T 3//2 which is consistent 
with a naive extension of our theory with the antiferro¬ 
magnetic fluctuations dominating [5] because of their in¬ 
creased effectiveness in momentum relaxation when com¬ 
pared to small q scattering. A detailed analysis is left for 
future work. Similarly a more detailed doping-dependent 
study of the Neel phase boundary Tjy(y) is necessary to 
compare with our predictions of a cross-over in power 
law. 

NbFe 2 is not unique in showing quantum multicritical¬ 
ity. YbRIi 2 Si 2 orders antiferromagnetically at low tem¬ 
peratures but the specific heat and Griineisen parameter 
at low temperatures obey power laws as would be ex¬ 
pected of a ferromagnet [2]. This could be the result of 
the presence of both ferro- and antiferromagnetic fluctu¬ 
ations [T9j[20]. However, this material is usually thought 
to he outside the Hertz-Millis scenario for quantum crit¬ 
icality because of Hondo breakdown effects Ell¬ 
in summary we have analysed the interplay of two 
quantum critical points in an itinerant magnet in both 
two and three dimensions. Our main prediction is that 
if quantum critical fluctuations of both ferro- and anti¬ 
ferromagnetic order are present then the specific heat, 
thermal expansion and Griineisen parameter will have 
the temperature dependence associated with just the fer¬ 
romagnetic modes. In addition, the correlation length 
of antiferromagnetic order will acquire the temperature 
dependence of the ferromagnetic correlation length, and 
this suppresses the boundary of the ordered phase (or 
region of applicability in two dimensions). We find that 
the boundary of the antiferromagnetic phase can undergo 
a crossover from its usual Hertz-Millis power law at low 
temperatures to the power law usually associated with 
a ferromagnetic instability at higher temperatures. We 
find no difference between a quantum bicritical and a 
quantum tetracritical point, as under renormalization the 
system always flows to the Gaussian fixed point where the 
interactions are zero. 
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